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RESUMEN 

El modelo estandar para la emision libre-libre de vientos ionizados se basa 
en la suposicion de un viento estacionario, isotropico y homogeneo. Sin embargo, 
durante la ultima decada estas suposiciones han sido cuestionadas. Existen objetos 
cuyo viento no se comporta como lo predice el modelo estandar. En este trabajo 
presentamos resultados para tres fuentes brillantes: P Cyg, Cyg OB2 No. 12 y WR 
147. Se ha reportado que estos objetos poseen vientos que se desvian de lo esperado. 
Reportamos el flujo observado y derivamos el tamano, el indice espectral y la tasa 
de perdida de masa para cada uno de ellos. Estos parametros nos permiten detectar 
la presencia de posibles asimetrias, inhomogeneidades y variaciones en la densidad 
de flujo. Estas desviaciones confirman un comportamiento distinto al clasico en 
estos vientos. 

ABSTRACT 

The classical model for free-free emission from ionized stellar winds is based 
on the assumption of a stationary, isotropic and homogeneous wind. However, since 
there exist objects whose wind behaviour deviates from the standard model, during 
the last decade these assumptions have been questioned. In this work, we present 
results for 3 bright sources: P Cyg, Cyg OB2 No. 12 and WR 147. These objects 
have been reported to possess winds that deviate from the basic assumptions. We 
have obtained flux densities, sizes, spectral indices and mass loss rates for each 
of the targets. These parameters allow us to analyze possible asymmetries, inho- 
mogeneities and time variations in the flux densities. These features confirm the 
nonclassical behaviour of these winds. 

Key Words: STARS: INDIVIDUAL (WR147, P CYGNI, CYGNUS OB2 
NO. 12) STARS: MASS LOSS 



1. INTRODUCTION 

Since winds from massive stars are ionized by the 
ultraviolet radiation from the underlying star, they 
can be studied observationally at radio wavelengths 
through their free-free emission. Until recently, mas- 
sive star winds were treated as isotropic and homoge- 
neous flows with constant velocity and known elec- 
tron temperature, degree of ionization and chemi- 
cal composition. In this way, based on the classic 
work of Panagia & Felli (1975) and Wright & Barlow 
(1975), it was possible to obtain a value for the mass 
loss rate M from radio observations. However, since 
the work of Abbott, Bieging & Churchwell (1981), 
who found that in some cases these assumptions are 
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not valid, various authors have found firm evidence 
for deviations from these basic assumptions. 

The existence of possible inhomogeneities has 
been studied, both observationally and theoretically. 
Radio observations seem to detect (marginally) 
blobs of material moving with the wind (Skinner 
et al. 1998; Contreras & Rodriguez 1999; Exter et 
al. 2002). In the optical region, peaks superposed 
on emission lines are suggested to represent individ- 
ual blobs (Moffat & Robert 1994). Theoretically, 
Cherepashchuck (1990) studied the discrepancy be- 
tween the expected and the observed X-ray luminos- 
ity. He concluded that a model including inhomo- 
geneities can explain this discrepancy and suggests 
that as much as 80% of the mass of the wind could be 
in the form of blobs. In fact, Williams et al. (1997) 
reported firm evidence of an asymmetric wind in WR 
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TABLE 1 
SOURCE SIZES AND FLUXES 



Source 



"0.7cm 

(PA) 
[arcsec] [deg] 



#3. 6cm (PA) 

[arcsec] [deg] 



#6cm(PA) 

[arcsec] [deg] 



So. 7 cm 

[mJy] 



$3. 6cm 

[mJy] 



~6cm 

mJy] 



P Cyg 0.06 x 0.03 (42) 0.19 x 0.13(160) 0.22 x 0.19(29) 9.7 ±2.1 8.0 ±0.2 6.0 ±0.1 

Cyg OB2 No. 12 0.04 x 0.03(32) 0.13 x 0.12(103) 0.24 x 0.12(176) 9.0 ±1.5 5.9 ±0.1 4.2 ±0.1 

WR 147 S 0.06 x 0.05(22) 0.25 xO. 21(120) 0.46 x0.35(163) 36.6 ±2.2 24.5 ±0.1 22.1 ±0.2 

WR147N ••• 0.37 x 0.30(110) 0.49 x 0.21(93) ••• 9.3 ±0.2 10.8 ±0.2 



Note.— Deconvolved size errors are ~ 0"01. 

147, making the WR star of this binary system a very 
good example of a non-classical wind source. 

Variable radio emission has been detected from 
several massive stars including P Cyg (Abbott, 
Bieging & Churchwell 1981; Contreras et al. 1996; 
Williams et al. 1997). This variability suggests that 
we are not dealing with a stationary wind. It has 
been proposed that these variations are due to a 
changing mass loss rate. However, the variability in 
the observed optical emission lines can be explained 
by the presence of an asymmetric stellar wind in ro- 
tation. In the case of a binary system, where there is 
an interaction surface between the two stellar winds, 
the variability can be due to the orbital sweeping of 
this interaction surface even in the absence of intrin- 
sic variability of the stellar winds (Girard & Willson 
1987; Luchrs 1997; Georgiev & Koenigsberger 2002). 

The study of the main assumptions of the classi- 
cal wind model is very important because mass loss 
rate values are usually derived assuming a classical 
thermal wind. Thus, in order to determine reliable 
mass loss rates, we should take into account possi- 
ble deviations (hacer algo en cada fucnte). In this 
paper, we present new VLA observations at 0.7, 3.5 
and 6 cm of three bright radio sources: P Cyg, Cyg 
OB2 No. 12 and WR 147. 

2. OBSERVATIONS 

We present two sets of observations taken with 
the Very Large Array (VLA) of the NRAO 3 on 1999 



3 The National Radio Astronomy Observatory is operated 
by Associated Universities Inc. under cooperative agreement 
with the National Science Foundation 



June 28 and September 3 at 0.7, 3.6 and 6-cm. Dur- 
ing both epochs the array was in the "A" configu- 
ration giving angular resolutions of ~ 0"05, ~ 0''2 
and ~ 0"4 for each wavelength, respectively. Dur- 
ing the first observing run all three sources were ob- 
served, while during the second only WR 147 was 
observed. The total on-source integration time was 
~ 1 hour for each source for the first run and 1.5 
hours for the second. Amplitude calibrators were 
1328±307 and 0137±331 and phase calibrators were 
2005±403 and 2015±371 on each observing run. For 
the first season, the observed bootstrapped flux den- 
sities for 2005+403 were: 1.18±0.06, 2.01±0.02 and 
2.54±0.03 mJy at 0.7, 3.6 and 6 cm respectively and 
for 2015+371: 1.35±0.09 at 6 cm. For the second 
run only 2005+403 was observed and a flux density 
of 1.99+0.01 at 3.6 cm was obtained. 

Data reduction was performed using the Astro- 
nomical Image Processing System (AIPS) software of 
the NRAO. At 3.5 and 6 cm we have followed stan- 
dard VLA procedures for editing, calibrating and 
imaging. At 0.7 cm, we have applied a small correc- 
tion to the data for atmospheric extinction using a 
zenithal optical depth of r — 0.1. In this way, we pre- 
sume to have minimized the expected tropospheric 
opacity dependence of the aperture efficiency. 

CLEANed maps were obtained using different 
data weighting schemes according to the source. 
For P Cyg and Cyg OB2 No. 12, an intermediate 
weight between Natural and Uniform was used (RO- 
BUST=0, Figures 1 and 2). While for WR 147, pure 
uniform weighting was preferred (Figure 3). In or- 
der to obtain flux densities and source sizes, we have 
used the AIPS IMFIT procedure. Based on a 2D- 
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Fig. 1. VLA maps of P Cyg and Cyg OB2 No. 12 at all observed wavelengths. All maps were obtained using AIPS 
routine IMAGR with an intermediate data weight (ROBUST=0). In the 3.6 cm map of P Cyg we can see a possible 
blob to the N-E in the wind, as well as one in the Cyg OB2 No. 12 map at 6 cm, to the S-E. Both weak emission 
regions are ~ 3 - 4cr above the map noise of 6.48 x 1(T 5 (P Cyg) and 5.72 x 1(T 5 (Cyg OB2 No. 12) at 3.6 and 6 cm, 
respectively. 



Gaussian fit to the source, we have determined inte- 
grated flux densities (Table 1) and deconvolved sizes 
(Table 2) for all sources. Distances and wind ter- 
minal velocity values were taken from the literature 
(Contreras et al. 1996); these data, together with 
the relation derived by Panagia & Felli (1975) and 
Wright & Barlow (1975), allowed us to derived mass- 
loss rates for each object (Table 2). Our values are 
not corrected for the effects of possible asymmetries 
or inhomogeneities, but they are expected to be cor- 
rect to an order of magnitude. 

3. RESULTS 
3.1. P Cyg 

This Bl la star has been studied in the radio since 
the early 1980's when Abbott, Bieging & Churchwell 
(1981) reported for the first time variable radio emis- 
sion. Radio emission variations are as large as ~50% 
without considerable changes in the spectral indices 
(Contreras et al. 1996). Several explanations related 
to wind parameters (ionization degree, wind veloc- 
ity, mass-loss rate) have been proposed, but none of 



them seems to fit the observed behaviour. In this 
work we have determined flux densities at 0.7, 3.6 
and 6 cm, based on our 2D Gaussian fits. Compar- 
ing these new flux densities with previous ones, we 
found that they have decreased by ~70%, ~35% and 
^34%, respectively, from their 1995 values (Contr- 

TABLE 2 
DERIVED PARAMETERS 



Source 



ao.7-3.6cm 



^3.6 — 6cm 



M .7cm 
[1CT 5 M ( 



©J 



P Cyg 0.12 ±0.13 0.50 ±0.05 0.5±0.2 

Cyg OB2 No. 12 0.26 ±0.10 0.63 ±0.06 2.7±1.0 

WR 147 S 0.25 ±0.04 0.18 ±0.02 1.0±0.4 

WR 147 N • • • -0.25 ±0.04 
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Fig. 2. VLA 
intermediate 
the WR star. 



CLEANed maps of WR 147 at 0.7 and 6 cm. Maps were obtained using AIPS routine IMAGR with an 
data weight (ROBUST=0). In the 6 cm map we can see two faint small emitting regions to the south of 



eras et al. 1996). This kind of behaviour has been 
observed by Skinner et al. (1998) in their long-term 
deep radio study. They found that P Cyg shows im- 
portant changes in its flux density where variations 
of 20-50% are frequently observed. 

Regarding spectral indices, Skinner et al. (1998) 
found two different values corresponding to two dis- 
tinct regions around P Cyg: the inner region (about 
10" in radius) presents a spectral index of 0.8 while 
the outer region of the nebula shows a spectral in- 
dex close to zero. In our case the spectral index 
behaviour seems to be the opposite: the 0.7-3.6 cm 
spectral index, coming from a region closer to the 
star, is almost flat (Table 2) while the 3.6-6 cm spec- 
tral index, coming from regions farther out from the 
star, is consistent with the 0.6 value for a thermal 
wind. The 0.7-3.6 cm spectral index has been al- 
tered by the decrease of the 0.7 cm flux density. As 
Skinner et al. (1998) proposed, since different wave- 
lengths sample different regions in the wind, each one 
with a completely different time-scale variation, one 
can expect that flux density changes may be poorly 
correlated or entirely uncorrelated. Then, we are ob- 
serving a rapid change in the wind emission at 0.7 
cm ocurring in the inner region but not observed in 
the outer regions of the wind. Wind variability has 
been observed at other wavelengths. Recently, long- 
term optical photometric studies by van Genderen, 
Sterken & De Groot (2002) and de Jager (2000) have 
described various types of stellar oscillations that P 
Cyg undergoes which could be related to changes in 
the wind parameters. Thus, in light of these new 
studies, it would be useful to examine again the pro- 
posed explanations mentioned above. 



Regarding morphology, we present 0.7, 3.6 and 
6 cm maps shown in Figure 1. From these maps, 
P Cyg seems to have a non-isotropic wind at all 
three wavelengths being more obviously elongated 
at 0.7 cm where its major /minor axis ratio is 2 (Ta- 
ble 2). The observed wind asymmetry is consistent 
with the results of Skinner et al. (1998), who found a 
slightly asymmetric structure surrounding the com- 
pact stellar wind source. Moreover, they found that 
the nebula around P Cyg is very inhomogeneous or 
filamentary at various scales: a compact core of a 
few tenths of an arcsec in diameter, an inner nebula 
of ~ 16 — 18" in diameter and an outer nebula ex- 
tending to 90". Although, in our 0.7 and 6 cm maps 
we can hardly see what could be emission from a 
blob denser than the surrounding wind, evidence of 
a blob is more obvious in our 3.6 cm map. Here we 
can see emission at a 4cr-level of a small condensation 
to the NE of P Cyg. This inhomogencity is about 
0"6 away from the central star which indicates that 
at an intermediate scale between the central core and 
the inner nebula (as defined by Skinner et al.) the 
wind is also inhomogeneous. 

The derived mass-loss rate for P Cyg, while 
somewhat lower than that reported by Contreras et 
al. (1996), is still consistent within the error, and 
consistent with the value of 1 - 3 x 10~ 5 M Q yr _1 
reported in the literature (Skinner et al. 1998). Since 
our two values were obtained using the 0.7 cm flux 
density, to avoid possible non-thermal emission con- 
tamination (which showed a very important decrease 
in 1999), it is not surprising that the 1999 mass- loss 
rate is lower than the previous one. P Cyg's wind 
and associated nebula, require more detailed study 
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because the radio variability has yet to be explained. 

3.2. Cyg OB2 No. 12 

Although no variable radio emission has been re- 
ported for this source, our new VLA observations 
show some variation in its 0.7, 3.6 and 6 cm flux den- 
sities. Both, 0.7 and 3.6 cm fluxes have decreased by 
^60% and ~18%, respectively, while its 6 cm flux 
has increased by ~12%. Our flux densities were ob- 
tained from a 2D Gaussian fit to the source. These 
fluxes are shown in Tabic 1. The flux variations have 
caused the spectral indices to change as well. The 
3.6-6 cm index has decreased from those reported by 
Contreras et al. (1996), becoming more consistent 
with the classical value for a thermal wind of 0.6, 
while the 0.7-3.6 cm index has decreased to a rather 
flat value of 0.3±0.1. 

Maps at the three observing wavelengths are 
shown in Figures 1 and 2. Based on its morphol- 
ogy, Cyg OB2 No. 12 does not show any deviation 
from a spherical wind at 0.7 and 3.6 cm (Figures I 
and 2), however, at 6 cm its shape is obviously elon- 
gated (Figure I). From a 2D Gaussian fit, we have 
obtained its angular size at all three wavelengths; it 
is clear that at 6 cm, the source has a major to mi- 
nor axis ratio of 2 (Table 1). The contour map shows 
that the 6 cm emission does not have a smooth oval 
shape but rather shows a kind of enhaced emission 
region to the south. Furthermore, if the source had 
an asymmetric wind, its shape would be elongated at 
all wavelengths, which is not the case. Since free-free 
emission depends on density as S oc n^, we propose 
that this protuberance could be explained as a blob 
denser than its surroundings. Thus, the wind of Cyg 
OB2 No. 12 deviates from the standard assumption 
of a homogeneous wind. 

3.3. WR 147 

This very interesting system is now considered 
the archetype of colliding wind binaries. It clearly 
shows the presence of a wind interaction region curv- 
ing away from the Wolf-Rayet(WR) star, and emit- 
ting non-thermal radiation (north radio component). 
Its southern radio component (WR star) shows time- 
variable radio emission as well as inhomogencitics 
(Contreras & Rodriguez 1999; Watson et al. 2002). 
In this work, we present new maps obtained from our 
0.7, 3.6 and 6 cm data (Figures 2 and 3). In these 
maps we can see both components and some of their 
fine structure. 

Flux density variability has been reported for 
both radio components (Williams et al. 1997; Se- 
tia Gunawan et al. 2001). Although Setia Gunawan 



et al. (2001) find variability at 6 and 21 cm for the 
north component of the system, from our new values 
(Table 1), we cannot see any significant variation in 
the 3.6 cm flux density of WR 147N, compared to 
its 1996 values. Thus, its 3.6 cm flux density may 
have remained constant, within the error, over an 
interval of ~3 yr. On the other hand, the WR 147S 
3.6 cm flux density shows a decrease of ~14% from 
its 1996 value. Since the 3.6 cm flux reported here 
is the combination of two observing runs for WR 
147, we can take this value as a highly reliable one. 
Then, this variation confirms the suggestion of WR 
147S variability proposed by Williams et al. (1997) 
and the variation reported previously by Contreras 
& Rodriguez (1999). However, this latter paper re- 
ports an increase of the radio emission while the 
present work reports a decrease in the same radio 
flux density. As Watson et al. (2002) have suggested 
for the variations in the radio flux at 6 cm, the pres- 
ence of inhomogeneities in the wind of the WR star 
can cause quite random variations in the radio emis- 
sion. A precise comparison of our 6 cm flux density 
with the 1995 value (Contreras et al. 1996) cannot 
be made because those data were taken at a lower 
angular resolution. However, this new value is con- 
sistent with the 1998 flux reported by Watson et al. 
(2002) at the same wavelength. Thus, it may be that 
flux variations observed at 3.6 cm, which traces ma- 
terial inner in the wind, cannot yet be observed at 6 
cm, which traces the outer regions. 

Morphologically, the high resolution MERLIN 
observations of Williams et al. (1997) and Watson 
et al. (2002) showed that at 6 cm both sources, WR 
147S and WR 147N, are asymmetric. However, Con- 
treras & Rodriguez (1999) found no clear elongation 
in the wind of WR 147S while the northern source 
was clearly asymmetric. With our new VLA ob- 
servations we still do not find any clear asymmetry 
in the southern source, and again the shocked re- 
gion (northern component) is obviously elongated. 
As discussed by Contreras & Rodriguez (1999), al- 
though we confirm that the northern source is elon- 
gated in the E-W direction, the dimensions derived 
using MERLIN data are smaller than ours because 
they resolve out extended emission. In Table 1 we 
report the deconvolved sizes for both sources at all 
three observed wavelengths. At present, it is quite 
clear that both radio sources show inhomogeneities. 
Based on MERLIN observations at 6 cm, Watson et 
al. (2002) have found that both radio components 
show radio emission structure. In fact, both emis- 
sion structures are highly variable as in the case of P 
Cyg's wind, but on a larger time-scale. In our maps 
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Fig. 3. VLA CLEANed map of WR 147 at 3.6 cm. The map was obtained in the same way as those of Figue 2 with 
a data weight of ROBUST=-l. We can see two blobs of emission to the south of the WR star and a small peak of 
emission to the E part of the shock region, see text. We show a set of fitting curves at different inclination angles. The 
best fit curve is that at 30°. 



we can see a slightly independent emitting region at 
the N-E side of the northern source that has not been 
detected in previous VLA observations. Giving the 
high resolution structure observed with MERLIN, 
we suggest that this small region could be a clump 
entering the colliding-wind zone as suggested by Se- 
tia Gunawan et al. (2001). Besides, we detect, at a 
3-4(7 level, two quite conspicuous condensations in 
the southern radio component located to the South 
at 3.6 cm (Figure 3). Besides, the clump directly 
to the south seems to be have a counterpart at 6 
cm (Figure 2), and additionally in the 6 cm map we 
can see another clump in the S-W direction, which 
does not coincide well with one any of the two de- 
tected at 3.6 cm. What is most surprising about 
the 3.6 cm blobs is that they seem to lie almost in 
the same direction and position as those reported by 
Contreras & Rodriguez (1999), based on data taken 
3 years earlier. Since we do not expect them to be 
the same ones previously observed, then it may be 
that as Watson et al. (2002) have suggested, the WR 
wind is directed in a north-south preferential outflow 
axis. 



Derived spectral indices are shown in Table 2. 
The northern source possesses a negative index in- 
dicative of non-thermal emission, which is consistent 
with previous results. Based on these values, we can 
see that both spectral indices for the southern com- 
ponent are lower than the expected value for a classi- 
cal thermal wind. This behaviour can be explained if 
we consider the extreme flux variability at all wave- 
lengths. 

3.4. Fitting the WR U7 Bow Shock 

We have previously fitted a colliding wind model 
to the radio morphology of WR 147N (Contreras 
& Rodriguez 1999). Given the improved quality of 
our new maps, we now rederive parameters for the 
system. In order to apply an analytic model, we 
will neglect inhomogeneities and possible mass loss 
variations, as well as the orbital motion of the driv- 
ing stars. The latter assumption is valid because 
the wind speeds are much greater than the orbital 
speeds. A thin shell model for the wind collision has 
two fundamental parameters under the assumption 
of steady, isotropic winds: the ratio of momentum 
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loss rates 77, and the angle <f> between the symme- 
try axis of the collision surface and the plane of the 
sky. We adopt the exact solutions of Canto, Raga & 
Wilkin (1996) for the collision surface, and apply the 
projection method described by Wilkin (1997) (see 
also, Raga et al. 1997). Along the line joining the 
centers of the two stars, the winds collide head-on 
at the stagnation point. If the stars are separated 
by a distance D, and the stagnation point is at dis- 
tance R from the source of the weaker wind, 77 may 
be obtained from the ratio R a /D = -q 1 / 2 / (1 + 77 1 / 2 ). 
This ratio is independent of the angle at which we 
view the system. In our fit, we assume that the peak 
emission of WR 147N corresponds to the stagnation 
point, and measure its position relative to WR 147S, 
0"556 ± 0"020 The observed separation of the two 
stars, 0"635 ± 0!'020 was taken from Williams ct 
al. (1997). We obtain r\ = 0.02 ± 0.01. This new 
value of 77 differs from the one obtained by Contr- 
eras & Rodriguez (1999) of 0.01, but is still consis- 
tent within error. Our revised value is also in agree- 
ment with that used by Dougherty et al. (2003) to 
reproduce the observed 5 GHz radio synchrotron and 
thermal emission. Our assumption that the bright- 
est part of the wind collision surface corresponds 
to the stagnation point is intuitive, since the nor- 
mal shock velocity is greatest at this location, but 
is also supported by simulated images of colliding 
winds (Raga et al. 1997). Given a value of 77, we 
then calculated the projected shock surface for sev- 
eral values of (j> (Fig. 3). Reasonable fits to the 
bow shock wings are obtained for 15° < (f> < 40°, 
with a best fit of = 30°. Thus, a lower limit 
to the inclination angle of the orbit would be 30°. 
Although our <f) = 30° curve is in good agreement 
with the northwestern wing of the bow shock, the 
northeastern wing appears to be more open, suggest- 
ing a non-axisymmetric morphology. Although ana- 
lytic models may be generated for non-axisymmetric 
bow shocks (Wilkin 2000), such models introduce 
too many free parameters to attempt a unique so- 
lution given the present maps (see Gaensler 2002). 
However, it would be interesting to constrain such a 
model with additional kinematic information derived 
from observations of spatially resolved line emission. 

4. SUMMARY 
All of our observed objects arc clear examples 
of non-classical winds. P Cyg deviates in at least 
two ways from the standard model: its wind is 
anisotropic and inhomogencous. Additionally, it is a 
highly variable source, both in its radio flux density 
at all wavelengths observed and in its morphologi- 
cal structure. This variability has been previously 



reported based on high angular resolution observa- 
tions. We confirm that at a lower resolution we still 
detect variations. We found that Cyg OB2 No. 12 is 
a non-spherical source at 6 cm and we suggest that 
the observed elongation is due to a blob denser than 
the surrounding wind material. Thus, we report Cyg 
OB2 No. 12 as an inhomogeneous wind source for the 
first time. WR 147 is one of the most interesting bi- 
nary systems. A highly structured wind has been 
reported at high angular resolution. In this work, 
we detect two possible clumps in the wind of the 
WR star. These blobs are almost in the same di- 
rection and position as those reported previously by 
Contreras & Rodriguez (1999). This confirms the 
suggestion of Watson et al. (2002) of a north-south 
preferential outflow axis. The shock region (north 
radio component) seems to possess a faint emission 
peak located to its N-E side which could be a clump 
crossing this zone. We have fitted a theoretical bow 
shock curve to the northern radio component based 
on our new VLA data. These data give a higher 
wind momentum ratio, 77, than the one derived pre- 
viously but it is still consistent within error. Our new 
best fitting curve corresponds to an angle (f> of 301*5 
degrees which gives a lower limit to the inclination 
angle of the orbit. 

In this way, our new VLA data confirm the non- 
standard wind behaviour of these three interesting 
sources. 
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